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The linear muffin-tin-orbital method combined with density lunctional theory (local approximation) 
and the semiempirical method of the intermediate neglect of the differential overlap (INDO) based 
on the Hartree-Fock formalism are used for the study of the F centers (O vacancy with two electrons) 
in cubic and orthorhombic ferroelectric KNbOa crystals. Calculations for 39-atom supercells show 
that the two electrons are considerably delocalized even in the ground state of the defect. Their 
wave functions extend over the two Nb atoms closest to the O vacancy and over other nearby atoms. 
Thus, the F center in KNbOa resembles electron defects in the partially-covalent Si02 crystal (the 
so-called E'l center) rather than usual F centers in ionic crystals like MgO and alkali halides. This 
covalency is confirmed by the analysis of the electronic density distribution. Absorption energies 
were calculated by means of the INDO method using the A self-consistent-field scheme after a 
relaxation of atoms surrounding the F center. For the orthorhombic phase three absorption bands 
are calculated to lie at 2.72 eV, 3.04 eV, and 3.11 eV. The first one is close to that observed under 
electron irradiation. For the cubic phase, stable at high temperatures, above 708 K, only the two 
bands, at 2.73 eV and 2.97 eV, are expected. 

PACS numbers; 77.84.Dy, 71.15.Fv, 71.10.-w, 77.80. Bh 



I. INTRODUCTION 

Due to its technological importance potassium nio- 
bate, KNbOs, a perovskite-type ferroelectric material, 
has lately been subject to numerous ab initio electronic 
structure calculations. In the cubic KNbOa crystal each 
O atom is surrounded by four K atoms, two Nb atoms 
and eight next-nearest O atoms. Many of the calculations 
were based on tlie local density approximation (LDA) 
combined eitherEJ'0 with the linearized muffin-tin orbital 
(LMTO)u or with the pseudopotential methodaa as well 
as with the linearized augmented plane waveETEl (LAPW) 
scheme. 

Complementary to this approach is the Hartree-Fock 
(HF) formalism. Compared to the LDA, the HF scheme 
has the advantage of the exact treatment of exchange 
interactions. Recent implementations have no restric- 
tions on the spatial form of the potential, no poten- 
tial effects due to use of mufhn-tin boundary conditions 
and/or space-packing empty spheres. It gives the effec- 
tive charges and suggests a bond-population analysis be- 
tween pairs of atoms, and last, it allows one easily to 
perform the calculation of excited states and optical ab- 
sorption energies. 



Since such calculations are quite time consuming, there 
exist only a few HF studies for perovskite systems; see, 
for example, the cluster calculations of Rcf. ^ Instead, 
a simplified (semiempirical) version of the HF method 
widely |kH|ewn as intermediate neglect of the differential 
overlapt2ltll (INDO) has been applied successfully tp-cal- 
culations for many [Oxide crystals—including MgO ,113 a- 
AI2O3 (corundum) ,113 and Ti02 .r % Tn recent studies of 
pure KNbOa and KTaOs crystalslijilj their electronic 
structure and equilibrium ground-state structure for sev- 
eral ferroelectric phases as well as F-phonon frequen- 
cies were reproduced in surprisingly good agreement with 
both LDA calculations and available experimental data. 

It is well understood now that point defects play an 
important role in the electro-optic and nonlinear r^ti- 
cal applications of KNbOa and related materials.llll In 
particular, reduced KNbOs crystals containing oxygen 
vacancies reveal short-pulse excitations which could be 
used for developing fast optical correlators Jl3 Its use for 
light frequency doubling is seriously affected by pres- 
ence of unidr;ntified defects responsible for induced IR 
absorption.Ej The photorefractive effect, important in 
particular for holographic storage, is also well known to 
depend on the presence of impurities and defects. 



1 



One of the most common defects in oxide crystals is 
the so-called £ center, an O vacancy Vq which traps 
two electronsBJ In electron-irradiated KNbOa a broad 
absorption band is observed around 2.7 eV atr*oom tem- 
perature and tentatively ascribed to F center£l (see also 
Rcf. [2^ ). This defect is also of theoretical interest for two 
reasons. 

(i) Due to a low local symmetry of the O sites in 
the lattice, the threefold-degenerate 2p-type excited state 
could be split into several levels responsible for several 
absorption bands. This effect has been observed a long 
time ago for the F'^ punters in corundum, but theoreti- 
cally it was examinecS only very recently. Upon cooling 
from a high temperature, KNbOa undergoes a sequence 
of phase transitions from a paraelectric cubic phase to 
ferroelectric tetragonal phase (at 708 K), then to the or- 
thorhombic structure (at 498 K), and finally to the rhom- 
bohedral (at 263 K) phase. The atomic positionajp all 
these phases have been determined experimentally.Ej Un- 
der these phase transitions the local symmetry of the O 
vacancy also changes, which can, in principle, affect the 
optical properties of the F centers. This problem has 
never been addressed. 

(ii) Qualitative theoretical analyses of the F centers in 
perovskites predict the effect of the symmetry breaking 
of one-electron orbitals associated with the asymmetric 
electron density delocalization over the_.two Nb atoms 
closest to the O vacancy: Nbi-yo-Nb2.L3 

To answer these questions, as well as to check the as- 
signment of the 2.7 eV absorption band, we study in the 
present paper the F center in KNbOa using the super- 
cell model and two different theoretical techniques: full- 
potential LMTO and INDO. 

The paper is organized as follows. In Sec. || we de- 
scribe the essential features of the methods used. In Sec. 
Ill we discuss the main findings obtained for the ground- 
state properties of the F center while Sec. IV describes 
the absorption energies calculated by means of the INDO 
method for two different (the cubic and orthorhombic) 
phases. Last, the main results are summarized in Sec. M 
where results obtained by the two (quite different) meth- 
ods are compared, and the effect of chemical bonding 
covalency on the F center properties is discussed. 



II. METHODS USED 

A. Local density approximation 

Even the crudest approximation, the LDA, to the den- 
sity functional theory has been successfully applied to 
predict structural and dynamical properties of a large 
variety of materials. Equilibrium volumes, elastic con- 
stants, phonon frequencies, surface reconstruction, and 
magnetism are just some examples of properties which 
could be successfully calculated for systems without par- 
ticularly strong electron correlations within the LDA (or 



LSDA, the local spin-density approximation). The LDA 
usually leads to some overbinding in solids (equilibrium 
volumes are typically 1-3 % underestimated). 

The LDA calculations for KNbOs performed in Rcf. 
^ yield an equilibrium volume which is « 5% too small, 
indicating that the LDA overbinding in this case is not 
considerably exceeding the "acceptable" limits. This is 
why in the present paper, although one may expect that 
a more accurate treatment of coixelation effects may be 
needed for some ionic compounds]^ we apply the LDA to 
the F center in KNbOs. The LDA exchange-correlation 
contribution is accounted_for by means of Perdew and 
Zunger's parasietrizationEZl of the calculations by Ceper- 
ley and Alder .c3 

The self-consistent solution of the one-electron equa- 
tion is performed by means of the LMTO method.H We 
have used the "atomic-spheres-approxipiation"|j (ASA) 
as weh as a "full-potential formalism"N (LMTO-FP). 
Whereas the LMTO-ASA uses potentials and charge 
densities that are made spherically symmetric inside 
(slightly) overlapping atomic spheres, no shape approx- 
imations are made in the LMTO-FP. The atomic re- 
laxations around the F center cannot be calculated by 
means of the ASA. We therefore performed the structural 
optimization by minimizing the LMTO-FP total energy 
calculated for a supercell. A similar method was used in 
earlier LMIO-FP simulations of defects in KCl (Ref. |3^) 
and MgO.Eil The supercell used in the present work con- 
tains 40 atoms for the perfect KNbOs (8 formula units) 
and 39 atoms plus one empty sphere in the F center case. 



B. INDO 

The INDO calculation scheme and the computer code 
CKUSTERD were discussed in detail in Rcfs. |o|jl| and 
p^ . With this code it is possible to perform both cluster 
and periodic system calculations containing hundreds of 
atoms as well as to carry out automated geometry op- 
timization which is especially important in defect cal- 
culations. In the periodic calculations the so-called large 
unit cell (LUC) model is used.E3 Its idea is to perform the 
electronic structure calculations for an extended unit cell 
at the wave vector k=0 in the narrowed Brillouin zone 
(BZ) which is equivalent to band calculations at several 
special points of the normal BZ, transforming to the nar- 
row BZ center after the corresponding extension of the 
primitive unit cell. In the KNbOs case the unit cell con- 
tains five atoms whereas the 2x2x2 extended (super)cell 
consists of 40 atoms. A detailed analysis of the KNbOa 
parametrization for the INDO method is presented in 
Ref. In that work considerable covalency was found 
of the chemical bonding in pure KNbOs. The effective 
charges found from MuUiken population analysis are (in 
units of |e|) -hO.543 for K, -f2.019 for Nb, and -0.854 
for O, which is very different from the expectation of 
the generally accepted ionic model: -1-1, -1-5, and —2, re- 
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spectively. This is in agreement with the effective atomic 
charges found in an experimental study of LiNbOa-Ej Our 
results emphasize a high degree of covalency of the Nb-0 
bond as may be expected from intuitive electronegativity 
considerations and the fact of a strong overlap between 
O 2p and Nb Ad orbitals and partial densities of states. 
We discuss below how covalency of the chemical bond- 
ing may have important consequences for the physics of 
defects in ferroelectrics, in particular for the F centers. 

To simulate F centers, we started with a 40-atom su- 
percell with one of the O atoms removed. In the cubic 
phase all O atoms are equivalent and have local sym- 
metry Ciy whereas in the orthorhombic phase there are 
two kinds of nonequivalent O atoms whose symmetry is 
lower, C2v or Cs (see below). After the O atom is re- 
moved, the atomic configuration of surrounding atoms is 
reoptimized via a search of the total energy minimum as a 
function of the atomic displacements from regular lattice 
sites. Calculation of the adiabatic energy curves for the 
ground and excited states permits us to find the optical 
absorption energy using the so-called A self-consistent- 
field (SCF) procedure according to which the Eats sought 
for is the difference of the total energies for the ground 
and excited state with the defect geometry of the ground 
state unchanged (vertical optical transition). 

To extend the basis set in the F center calculation, ad- 
ditional Is, 2p atomic orbitals were centered on the O va- 
cancy. Their parameters were chosen close to those used 
in the F center calculations in MgO crystaha the or- 
bital exponent C(ls)=0.65 (a.u.^^), C(2p)=0.50 (a.u"!), 
the relevant electronegativities are zero (Is) and —3 eV 
(2p), and the bonding parameters /3=0 for both Is and 
2p orbitals, respectively. During the defect geometry op- 
timization, we make no a priori assumptions on the elec- 
tron density distribution. 



III. GROUND-STATE PROPERTIES 
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FIG. 1. (a) Partial F-center density of states for KNbOa 
calculated within the LDA using a 39-atom supercell. Ener- 
gies are given relative to the Fermi level. This band structure 
is metallic because the finite size of the supercell causes the 
defect band to have a width of ~0.8 eV so that it overlaps 
with the too- low- lying conduction band ("LDA gap error"), 
(b) As in (a) but with application of an external potential that 
upshifts the Nb d states in energy by «2.5 eV. Zero energy 
corresponds to the top of the F center band. 



A. LDA calculations 

The band structure as derived from the straight LDA 
underestimates the gap between occupied and empty 
states. Since the supercell which we use is rather small 
(40 atomic sites) , the defect states form a band of a finite 
width («0.8 eV). These two effects cause the defect band 
to overlap with the conduction band, and the supercell 
calculations within the LDA predict KNbOa with the F 
centers to be a metal. This is seen from Fig. 11(a). This 
affects the charge distribution. The number of electrons 
in the vacancy "atomic sphere" is 0.24, much smaller 
than obtained from the INDO (0.6; see below). We can 
artificially increase the gap by applying an upshift in each 
iteration to the Nb d states. Figure H(b) shows such a cal- 
culation, and now the defect band lies completely within 
the gap. The energy scale was chosen so that the top of 
the defect band (the highest occupied state) is at £'=0. 



The band plot. Fig. || also shows that the defect band 
disperses over «0.8 eV due to the small supercell size. 
After adjustment of the gap the (self-consistent) calcu- 
lation yields 0.6 electrons inside the O vacancy sphere, 
i.e., more than twice the amount found before and very 
close to the INDO calculation for the relaxed structure. 
When relaxations (see below) are included, the LMTO 
calculation yields a lower electron number in the vacancy 
sphere. This is simply caused by the outward motion of 
the nearest neighbors (Nb). 

The relaxation of atoms surrounding the F center was 
first calculated within the LDA without any attempt to 
correct for the effect of the overlap between the defect 
state band and the conduction band. First, the nearest- 
neighbor Nb atoms were relaxed, and the result is il- 
lustrated in Fig. ||(a) which shows the total energy as 
a function of the outward displacement, A^, of the Nb 
atom from its equilibrium position in the bulk crystal. 
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TABLE I. Coordinates of three kinds of 14 atoms surrounding F center placed at the coordinate origin in the cubic phase 
(in units of the lattice constant ao=4.016 A). 



Atom 




Lattice coordinates 






Displacement 






without additional AO 






Nb 








±(| + A.) 




= 0.08320 


K 


±(1 + A,) 







A^ 


= 0.01324, Ay = 0.01324 





±(5 - A.) 





±(i-A.) 


A^ 


= 0.02372, A^ = 0.00447 


O 





±(5 -A.) 


±(5 -A.) 


A« 


= 0.02372, A^ = 0.00447 






with additional AO 






Nb 








±(5 + A.) 




A;, = 0.06550 


K 


±(i+A.) 







A^ 


= 0.00920, Aj; = 0.00920 


O 


±(5 -A.) 





±(5 -A.) 
±(|-A.) 




= 0.01920, A^ = 0.00334 








±(^-A.) 


A« 


= 0.01920, A^ = 0.00334 




FIG. 2. Calculated band structure for the 39-atom super- 
cell modeling the F center. The Nb d states are shifted up as 
in Fig. 1(b). The horizontal line {"Ef") indicates the top of 
the defect band, the highest occupied state. 

The relaxed atomic positions correspond to A2=3.5% of 
the lattice constant of ao=4.016 A. This is about half 
the relaxation found in the INDO calculation (see be- 
low). Further, the relaxation energy found here, « 0.5 
eV, is much smaller than the value of 3.7 eV obtained in 
the INDO calculation. 

We do not wish to rely on total energy calculations 
where we applied 2.5 eV upshift to the Nb d bands. How- 
ever, self-consistent calculations using supercells large 
enough for obtaining a small width of the defect band 
nonoverlapping with the LDA conduction band are im- 
practical. An approximate calculation was instead made 
of using the same supercell size as before but sampling 
only the F point of the BZ in the k-space integration. 
This is the point where the defect band has its mini- 
mum energy and lies inside the gap, even in the LDA 
calculation. Further, this sampling is the same as used 
in the INDO calculation described below [the so-called 
large unit cell (LUC) model]. As expected, this changes 
the charge distribution and the atomic relaxations. The 
value of 4.8% [Fig. |l)(b)] is closer to that obtained 
in the INDO calculation. The relaxation energy 1.2 eV 




Az{Nb) (%ofa) 

FIG. 3. (a). LMTO-FP calculation within the LDA of the 
total energy vs the outward displacement, 5z, of the two Nb 
atoms closest to the O vacancy. Note that no corrections were 
made in this case to ensure that the defect band lies entirely 
inside the gap. Integration in k space used 40 points in the 
irreducible part of the BZ. (b). Same as (a) but only a single 
r point was here included in the k-space sampling, in order to 
simulate the LUC used in the INDO calculations. The curve 
is a high- order polynomial fit. 
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TABLE II. Atomic coordinates for a perfect orthorhombic phase of KNbOa [in terms of lattice parameters a=3.973 A , 
6=5.695 A , c=5.721 A (INDO calculated in Ref. as well as positions of the two Nb atoms nearest to the F center related 



to a 0(1) vacancy (point symmetry C2v) and 0(11) vacancy (symmetry Ca, respectively). The sketch of the orthorhombic unit 


cell is shown 


in Fie. W. 














Atom 


a 


b 


c 




A 












Perfect crystal 










K 








A. 


A, 


= 0.0209 






Nb 


2 





1 

2 










Oi 








i+A, 


A, 


= 0.0347 






On 


2 




i+A, A,=-0.0028 






A, 


=0.0347 


On 


2 




i+A, A,=-0.0028 






A, 


=0.0347 






F centers of C2v symmetry: no additional AO 










Nb(I) 







i+A, A.=0.08333 






A, 


=0.0132 


Nb(II) 


-i+A. 





i+A, A,=-0.08333 






A, 


=0.0132 






F centers of C2v symmetry: with additional AO 










Nb(I) 







i+A, A,=0.0656 






A 


— n 01 fifi 

— U.U -LUU 


Nb(II) 


-i+A. 





i+A, A,=-0.0656 






A, 


=0.0106 






F centers of Cs symmetry: no additional AO 










Nb(I) 


1 
^1 


Ay 


i+A, A,=-0.0468 






A, 


=0.02502 


Nb(II) 


2 


i+A, 


A, A,=0.0244 






A,= 


-0.03735 






F centers of Cs symmetry: with additional AO 










Nb(I) 


1 


Ay 


i+A, A,=-0.0361 






A, 


=0.0192 


Nb(II) 


\ 

2 


I+A, 


A, A,=0.0187 






A,= 


-0.0286 



is also somewhat closer to the INDO result. 



B. INDO 

1. Cubic phase 

The positions of 14 atoms surrounding the F center in 
a cubic phase after lattice relaxation to the miminum of 
the total energy are given in Table | for the two cases of 
with an additional atomic orbital (AO) on vacancy site, 
and without the AO. The conclusion is that the largest re- 
laxation is exhibited by the two nearest Nb atoms which 
are strongly (by 6.5% of ao) displaced outwards the O 
vacancy along the z axis. This is accompanied by a much 
smaller, 0.9% outward displacement of K atoms and by 
a 1.9% inward displacement of O atoms. The two Nb 
atoms give the largest (« 80%) contribution to the lat- 
tice relaxation energy (3.7 eV) whereas O atoms give the 
most of the rest energy gain of 1 eV. 

The analysis of the effective charges of atoms surround- 
ing the F center shows that of the two electrons associ- 
ated with the removed O atom only « — 0.6|e| is localized 
inside Vq and a similar amount of the electron density is 
localized on the two nearest Nb atoms. The F center 
produces a local energy level, which lies ~ 0.6 eV above 
the top of the valence band. Its molecular orbital con- 
tains primarily contribution from the atomic orbitals of 
the two nearest Nb atoms. 




FIG. 4. Orthorhombic phase of KNbOa: equilibrium ge- 
ometry. Atomic coordinates are given in Table O. 

2. Orthorhombic phase 

The orthorhombic phase of KNbOa is important since 
it is stable in a broad temperature range around room 
temperature and thus is subject to most studies and 
practical applications. The atomic positions for the 
perfect ortjaorhombic cell shown in Fig. ^ were calcu- 
lated earliertj and includedria Table |ll[ Their agreement 
with the experimental dataEJ is very good.ll3 

The displacements of Nb atoms nearest to Vo are given 
in the same Table || for the both kinds of F centers ex- 
isting in this phase and, again, with and without atomic 
orbitals centered at the O vacancy. In fact, Nb displace- 
ments are very similar in magnitude (6.6%) and also close 
to those found for the cubic phase. The relevant relax- 
ation energies are considerable (3.6 eV) and nearly the 
same as that for Nb relaxation found in the cubic phase. 
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TABLE III. Calculated absorption (iSaba, eV) and 
Nb-atom relaxation (-Brci, eV) energies for the F center for 
the cubic and the orthorhombic phases 



Symmetry, phase 








Era, 


C^v, cubic 


2.73 


2.97 




3.7 


Cs, orthorhombic 


2.56 


3.03 


3.10 


3.6 


C2v, orthorhombic 


2.72 


3.04 


3.11 


3.6 



IV. OPTICAL PROPERTIES 

Because of the C4v local symmetry of the F centers in 
the KNbOs cubic phase, its excited state splits into two 
levels, one of which remains twofold degenerate. Our 
ASCF calculations predict the two relevant absorption 
bands: at 2.73 eV and 2.97 eV (Table Neglect of 



the additional orbitals centered on the Vq slightly affect 
this result, the relevant absorption energies turn out to 
be 2.67 eV and 3.02 eV, respectively. 

Around room temperature in the orthorhombic phase 
there exist two kinds of F centers associated with two 
non-equivalent O atoms revealing the C2v and Cs sym- 
metry; cf. Table ^ The corresponding three absorption 



bands for each of them are shown in Table IIL Their 
difference is the largest for the lowest-energy bands (0.16 
eV) and very small for other two bands. 



We presented a strong argument that the 2.7 eV att. 
sorption band observed in electron irradiated crystalsEll 
could be due to the F centers, and predicted existence of 
two additional absorption bands (at 3.04 eV and 3.11 eV) 
for the same defect in the orthorhombic phase of KNbOa 
(see also discussion in Ref. |2]). At higher temperatures 
where the cubic phase is stable, the latter two energies 
which hardly could be separated experimentally because 
of the large half-width of absorption bands, degenerate 
into a single, double-degenerate level at 2.97 eV. 

Our results alsp-suggest that the blue-light-induced-IR 
absorption effectta mentioned in the Introduction could 
be triggered by the F center absorption which may lead 
to its subsequent ionization where an electron is trans- 
ferred to the conduction band. The retrapping of this 
electron by another defect is then responsible for the in- 
frared absorption when the retrapped electron is excited 
to the conduction band. We recall that the UV excitation 
(blue- light) energy used in these experiments is very close 
to our calculated absorption energy of the F centers. 

The reason for the discrepancy in LDA and INDO re- 
laxation energies for the F center needs further study; 
probably it arises due to strong electron correlation ef- 
fects which are crudely approximated in the two methods 
in quite different ways. 
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V. SUMMARY 

Our results are in a sharp contrast with what is known 
for F centers, in ionic crystals (in particular, in MgO and 
alkali halidesl3E3) where the two electrons are well local- 
ized by the Vq in the F center ground state. Evidently, 
this discrepancy arises from a considerable degree of co- 
valency of the chemical bonding in KNhO^ which is ne- 
glected in all previous models of defects in this material 
(as well as similar ABO3 perovskites, e.g. Refs. ^p6[ ); 
the only exception known to jis is an Xa cluster calcula- 
tion of F centers in LiNb03.E3 

Electron defects similar to what we have observed here 
are known, in particular, in partly covalent Si02 crystals 
(e.g., in the so-called E[ center an electron is also not 
localized inside Vq but its wave function mainly over- 
laps mith the sp^ orbital centered on the neighboring Si 
atonfa). 

We found that the ground state of the F center is asso- 
ciated with a strong symmetrical relaxation of the two 
nearest Nb atoms outwards relative to the O vacancy. 
These Nb atoms remain to be identical; i.e., we did not 
see formation of dipole moments of the Nbi-yo-Nb2 type, 
as suggested in Ref. Note that the relevant relaxation 
energy is several eV which is typical for many point de- 
fects in ionic and partly ionic solids. Its magnitude is by 
several orders of magnitude larger than the tiny energy 
gain due to the phase transitions (meV per cell). 
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